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We apply the dynamical systems tools to study the asymptotic properties of a cosmological model 
based on a non-linear modification of General Relativity in which the standard Einstein-Hilbert 
action is replaced by one of Dirac-Born-Infeld type. It is shown that the dynamics of this model is 
extremely rich: there are found equilibrium points in the phase space that can be associated with 
matter-dominated, matter-curvature scaling, de Sitter, and even phantom-like solutions. Depending 
on the value of the overall parameters the dynamics in phase space can show multi-attractor structure 
into the future (multiple future attractors may co-exist). This is a consequence of bifurcations in 
control parameter space, showing strong dependence of the model's dynamical properties on the 
free parameters. Contrary to what is expected from non-linear modifications of general relativity of 
this kind, removal of the initial spacetime singularity is not a generic feature of the corresponding 
cosmological model. Instead, the starting point of the cosmic dynamics - the past attractor in the 
phase space - is a state of infinitely large value of the Hubble rate squared, usually associated with 
the big bang singularity. 

PACS numbers: 04.20.-q, 04.50.Kd, 95.36,+x, 98.80.-k, 98.80.Bp, 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 



to be cured - is replaced by one of the DBI form: 



Attempts to modify the Einstein-Hilbert (EH) action 
of General Relativity (GR) 
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where R is the Ricci curvature scalar, and A- the cos- 
mological constant (k 2 — 8irnipf = 8ttG), have been 
motivated by a number of reasons. In particular, renor- 
malization at one-loop demands that the Einstein-Hilbert 
action be supplemented by higher order curvature terms 
l]. 1 Besides, when quantum corrections or string theory 
are taken into account, the effective low energy action for 
pure gravity admits higher order curvature invariants Q . 

There are additional ways to modify the EH GR ac- 
tion. For instance, the one based on the Dirac-Born- 
Infeld (DBI) procedure for smoothing out singularities 
d, [j|. 2 According to this procedure the original La- 
grangian density £ = \f\g~\L - whose singularities are 
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1 Higher order actions are indeed renormalizable but not unitary 
0- 

2 The proposal to remove initial as well as final singularities in 
modified gravity has been given in Ref. [6|]. It was shown there 



£ 



£ 



DBI 




where the scale /i sets the limit of energy density acces- 
sible to the theory. 

A combination of the above possible modifications, i. 
e., a DBI- type action containing an EH term plus ad- 
ditional higher curvature terms within the square root, 
could supply an additional cosmological scenario where 
to look for alternative explanations to several phenomena 
such as inflation and the present speedup of the cosmic 
expansion. 

Several theories of gravity of this kind have been pro- 
posed since long_ago in Q, and in more recent years, 



for instance, in [8( (see also @). To be phenomenolog- 
ically viable, these modifications have to satisfy several 
physically motivated requirements Q: 

1. Reduction to EH action at small curvature, 

2. Ghost freedom, 

3. Regularization of some singularities (as, for in- 
stance, the Coulomb-like Schwarzschild singular- 
ity), and 



that the addition of a R? term to otherwise divergent modified 
gravity makes it regular. 
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4. Supersymmetrizability. 

The later requirement is quite stringent and, for most 
purposes, might be excluded. A DBI-type theory that 
fulfills the above requirements can be based on the fol- 
lowing action Q: 

Sdbi = ^J d 4 x^/\g~\ (l - y/l - a K 2 R + /?K 4 g) , (1) 

where Q = R 2 - tR^W + R^wR^™ is the Gauss- 
Bonnet (GB) invariant. It has been demonstrated in 
Ref . Q that this action has the EH leading term at small 
curvature, besides, for an appropriate region in the pa- 
rameter space it shows indications for the cancellation 
of the Coulomb-like Schwarzschild singularity, and it is 
ghost-free. 

Given that such a regular theory of gravity could 
have impact on the cosmology of our universe, in this 
paper we aim at exploring the asymptotic properties 
of a Friedmann-Robertson- Walker (FRW) cosmological 
model based on such a non-linear modification of General 
Relativity. 3 Our study will rely on the use of the stan- 
dard tools of the dynamical systems (see, for instance, 
Ref. [Tl|). The present investigation represents a gener- 
alization of the study in Ref.[12| to include dependence 
not only on the Ricci scalar but also on the Gauss-Bonnet 
term, which is one of the few ways to overcome the pres- 
ence of ghosts. Notwithstanding, we will not make a 
comparison of the results we obtain in this paper with 
the results of [l2j|, since the latter refer to a completely 
different non-linear modification of general relativity hav- 
ing nothing in common with the DBI modification stud- 
ied here. 

The paper has been organized as it follows. In the 
next section the basic mathematics of the model (fTJ) are 
exposed, including the derivation of the relevant FRW 
cosmological equations. The main section III is devoted 
to the study of the asymptotic properties of the model. 
A detailed discussion of the results is given in section IV, 
while brief conclusions are given in the last section V. 
Here we use natural units (k 2 = 8ttG = h = c = 1). 

II. EQUATIONS OF MOTION 

We start by writing the action ([1} in a more compact 
form [|: 

Sdbi = ^ J d 4 x^W\ f(R,G), (2) 

where 



3 The asymptotic properties of a big class of f(R, Q) theories, with 
the aim to study future singularities has been discussed in [lOj . 



f(R, Q) = \ (l - V~l - ok?R + /3k 4 g) . (3) 

Written in the form @ the action (|T|) has the appear- 
ance of an f(R, (?)-theory of the kind studied, for instance 
in Ref. J 131], so that we can safely apply known results of 
these investigations. 

The equations of motion that can be derived from ([2]) 
- plus the addition of a matter action piece - by varying 
in respect to the metric g M „ are the following [13j : 

-\g^R) + 4(i? CT ,v CT v A1 / g + Rv^vjg 

-g^R* T V a V T fg) + 2[2R, iauT V a V T 
— RVpVvlfg + 2fg(RR^ l/ — 2R^ a R a v 
+2R aT R^otv + Rn<jTpRJ 7Tp ) = K 2 Tfj, u , (4) 

where f R = 5 R f, fg = dgf, and □ = 5 CTr V CT V r . Worth 
noticing that, in case the function / were a function of the 
curvature scalar only, i. e., in case / = f(R) (formally 
this case corresponds to the limit Q — > 0), the above 
field equations reduce to the equations of motion of an 
/(-R)-theory [3, [Hj]. For purpose of comparison with the 
equations of EH GR, equations ^ can be recast into the 
form that resembles Einstein's equations 13]: 

— Tflv-vR = K lff{T^v + T™ rv ), (5) 

where k 2 ^ = n 2 /(f R — 4Dfg) is the effective gravita- 
tional coupling, and we have introduced the following 
curvature (effective) stress-energy tensor acting as an ad- 
ditional source of the field equations: 

K 2 T curv = { y i y v _ g ^ U )f R + 1 S/U ,(/ - Rf R ) 

-2(2R^ T V°V T - RV^7 v )f g - 4(i? £n/ V CT V M / s 
+R^V a VJ g - ,g^i? CTr V ff V T /e) - 2fg(RR^ 
2Rp b aR G y + 2R aT Ry uGTV + R yuGT pRj TTp ). (6) 

The trace of equation (Q| leads to the following con- 
straint: 

3Df R + Rf R -2f + 2RDf g 

+2[G - 2R aT \7 (7 V T }fg = n 2 T. (7) 

For a spatially flat FRW universe metric 

ds 2 = -dt 2 + a{t) 2 (dr 2 + r 2 (dQ + sin 2 9dip 2 ), 

equations (jj) can be written in the form of the following 
set of cosmological equations: 
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6H 2 f R = Rf R -f- 6Hf R 

-24H 3 fg+gf g + Pm , (8) 
AHf R = 2Hf R + 2fg(12H 3 - G/3H) 

-2f R - 8H 2 f g - p m , (9) 
p m = -3Hp m , (10) 

where the dot accounts for derivative with respect to the 
cosmic time t, and for simplicity we are considering dust 
as background fluid (it could be, for instance, the dark 
matter component). 

Starting with equation ([3]) it is straightforward to ob- 
tain the following relationship between the derivatives of 
the function f(R, Q): 



y- = = A (11) 

JR JRR 

where the constant A = — K 2 (3/a — —p/a (recall that 
we use the units system where n 2 = 1). In the present 
paper, for dcfinitcness, we shall choose only non-negative 
a > and /3 > 0. In consequence, the constant A will be 
a negative magnitude. 

After (fTTj) . and recalling that 

G = 24H 2 (H + H 2 ), R = 6{H + 2H 2 ), 

we are able to rewrite equations (I8I9P in a more compact 
yet simpler form: 

6H 2 f R (l + 4 AH 2 ) = Rf R (l + 4 AH 2 ) 

-f -6Hf R (l + 4AH 2 )+p m , (12) 
4Hf R + 16AHHf R = 2Hf R (l + 4 AH 2 ) 

- Pm -2f R (l+4AH 2 ). (13) 

Recall that, according to ([5]), the equation (TT2"j) can be 
written in the form of a standard Friedmann equation: 



III. DYNAMICAL SYSTEMS STUDY 

The dynamical systems tools offer a very useful ap- 
proach to the study of the asymptotic properties of the 
cosmological models [lf|. In order to be able to apply 
these tools one has to (unavoidably) follow the steps enu- 
merated below. 

• First: to identify the phase space variables that 
allow writing the system of cosmological equations 
in the form of an autonomus system of ordinary 
differential equations (ODE). There can be several 
different possible choices, however, not all of them 
allow for the minimum possible dimensionality of 
the phase space. 

• Next: with the help of the chosen phase space vari- 
ables, to build an autonomous system of ODE out 
of the original system of cosmological equations. 

• Finally (some times a forgotten or under- 
appreciated step): to identify the phase space 
spanned by the chosen variables, that is relevant 
to the cosmological model under study. 

After this one is ready to apply the standard tools of 
the (linear) dynamical systems analysis (see, for instance 
Ref. H2|). 

We will split our study into two parts. First, in sub- 
section A, we explore the asymptotic properties of the 
simplest situation when in (fTJ) we set Q = 0. Although 
in this case the corresponding DBI theory has a ghost, 
spin-zero degree of freedom, so it does not admit a stable 
de Sitter background [Hj], this case is easier to handle 
and could signal to relevant asymptotic behavior with 
which to compare the results of the more general ghost- 
free case {Q ^0), in the second part of our investigation 
(subsection B). 

A. Einstein-Hilbert DBI Case (G = 0) 



3H - K e ff(p m + Pcurv), K eff - zf^l+^AH 2 )' 
tiffPeurv = f - 2fR{l l AAH2) ~ 3if-^. (14) 

Notice that written in this latter form, due to the fact 
that A < 0, the Friedmann equation blows up at H — 
±1/(2\^A) = ±y/a/fi/2 (it is in fact the effective grav- 
itational coupling K 2 e ff who does). However this is not 
pathological since the above is just a convenient way to 
write the original equations. 

Our goal will be to write the above cosmological equa- 
tions in the form of an autonomous system of ordinary 
differential equations, so that we could apply the stan- 
dard tools of the dynamical systems to find the equilib- 
rium configurations that could be associated with rele- 
vant cosmological solutions. 



1. Autonomous System of ODE 

Our aim here is to write the following system of cos- 
mological equations (basically equations (Tl2"|) . and (TIB")) 
with A = 0): 

6H 2 f R = Rf R -f- 6Hf R + Prn , 

4Hf R = 2Hf R - 2f R - Pm , (15) 

in the form of an autonomous system of ordinary differen- 
tial equation (ASODE). To that purpose let us, following 
Ref.[19j, to introduce the dimensionless variables x, y, z: 



4 Notice that our choice of phase space variables differs from that 
of Ref. |19(. in the definition of one of the variables (y). 



4 



TABLE I: Properties of the equilibrium points of the autonomous system (|17l) . 
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TABLE II: Eigenvalues of the jacobian matrices corresponding to the critical points in table [i] 
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-2 
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3 
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4 
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4 

-3 


6 
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3 
2 


-4 

%/2T 
2 



R 
6H 2 ' 



3aift 



Hf R 



(16) 



Starting from equations (|12|) , (fT3|) . (|16j) and, after a bit 
of algebra, it is uncomplicated to obtain the following 
autonomous system - the right-hand-side of the ODEs 
do not depend explicitly on the time parameter r -, that 
describes the cosmological dynamics of this model: 



the curvature source in the right-hand-side of Einstein's 
equations. 

The first equation in (|15[) can be written in the form 
of the following (Friedmann) constraint: 



1 = 



R 



f 



h 



6H 2 f R 6H 2 QJrH 2 Hf R ' 
or, after the above choice of phase space variables: 



x' = -2x(x - 2) - z(y - 2x) , 

y' = -2y(x-2), 

z = z 2 - xz + 2(x - 2) + 

+3(1 + x - y + \/y 2 - 2xy) . (17) 

Here the comma denotes derivative in respect to the time 
variable r = In a (basically the number of e-foldings of 
inflation). 

The following expressions will be useful: 



fl m = 1 + x — y — z + \/y 2 — 2xy = 1 — O c 



(18) 



1 — x, w 



eff 



where q = — 1 — H / H 2 is the deceleration parameter, 
and w e f f is the effective equation of state parameter for 



where we have conveniently defined Cl m = p m /6H fn, 
and flcurv is the dimensionless curvature (effective or 
"dark" ) energy density. 

The next step is to identify the phase space relevant for 
the cosmological model associated with the autonomous 
system of ODE (fl~7|) . Since we are considering positive 
a > =>• fii > 0, then Q m > 0. Besides, since the 
expression under the square root in the intermediate ex- 
pression in Q18J1 must be positive, so that f2 m be real, 
then the following conditions are relevant to define the 
phase space: 



s + vV - 2xy > , y>2x, y>0, (19) 
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where the condition y > 2x holds only for non-negative 
x > 0. For x < the expression under the square root is 
always positive. 

Our most precise definition of the corresponding phase 
space is the following: 



®A = {{x, y, z) : \jy 2 - 2xy >y + z-x-l, 

y>0, y>2x}.(20) 




FIG. 1: Probe paths in phase space originated from different 
sets of initial conditions for f(R) = (l — y/1 — an 2 R) . The 
flux in time is shown in the lower panels. The trajectories in 
emerge from the past attractor (decelerated expansion 
phase corresponding to the equilibrium point Pi = (0, 0, 1) in 
Tab. HJ, and are attracted into the future by the curvature 
dominated inflationary fixed point P4 = (2,4, —I). 



The past attractor in the phase space is the equilib- 
rium point Pi — (0, 0, 1) in Tab. Q] Although the corre- 
sponding cosmological phase is dominated by the curva- 
ture (flcurv — 1), the expansion is decelerated {q = 1). 
The effective curvature energy density mimics radiation 
(w e ff = 1/3). This state is characterized by infinitely 
large values of the Hubble parameter squared H 2 — > 00, 
and fa oc a(t) - the scale factor. 

The fixed point P 2 = (0,0,-1) in Tab. U is a sad- 
dle equilibrium point. It is associated also with decel- 
erated expansion (q = 1), but this time there exists a 
scaling between the energy density of the background 
fluid Q m — 2 and that of the effective curvature com- 
ponent Slcurv = — 1: Qm/Qcurv = — 2. The fact that 
the dimensionless energy density parameter for the dust 
f2 m > 1, should not bother us. Recall that we have cho- 
sen a convenient definition Q, m = p m /6H 2 /r, so that, 
provided that 2 fa < Q m = p m /3H 2 , we will be faced 
with f2 TO > 1. During this phase of the cosmic evolution 
the effective energy density of the universe also mimics 
radiation. 

The equilibrium point P3 = (—1,0, —3) is also a saddle 
and is associated with a scaling of both energy density 
components: £l m /Qcurv = —3/2. In this case the dynam- 
ics of the expansion is "super" -decelerated (q = 2), while 
the effective equation of state parameter mimics that of 
stiff fluid (w e ff = 1). 

The fixed point P4 = (2, 4, —1) is the future (late-time) 
attractor in ^ a- It is associated with a curvature domi- 
nated (Clcurv — 1) &m = 0), inflationary (q = —1) solu- 
tion. The effective curvature fluid mimics a cosmological 
constant A (w e ff = — 1). Another curvature-dominated, 
inflationary, with A mimicry, twin-solution is associated 
with the spiral saddle equilibrium point P§ = (2,9/2,0) 
in Tab. Q] 

The saddle fixed point P5 = (2,4,3) does not re- 
ally correspond to a meaningful cosmological scenario. 
Actually, in this case fi m = p m /6i? 2 /^ = —4, but 
since we are considering only non-negative a > 0, then 
Jr. = a/yjl — aR > 0. Therefore, even if Cl m is an effec- 
tive (convenient) parametrization, it can not be negative. 



2. Equilibrium points 

The next step is to find the roots of the system of equa- 
tions x' = , y' = , z' = 0. Then one linearizes the 
ASODE (fl~7|) - means that we expand (fT7|) in the neigh- 
borhood of the equilibrium point keeping only the linear 
terms - and finds the eigenvalues of the linearization ma- 
trix. 

The most relevant properties of the equilibrium points 
of the autonomous system of ODE (jTTJ) are summarized 
in table |H Table |TT] displays the eigenvalues of the respec- 
tive linearization or Jacobian matrices. 



B. Ghost-free Case (Q / 0) 

1. Autonomous System of ODE 

Now we turn to the more general case when in equa- 
tions (O, flTBl . the constant A ^ 0. In order to trans- 
form these equations into an autonomous system of ODE 
we introduce the following phase space variables: 



R 
6H 2 ' 



1 



3aH 2 
1 



Hfn 



AAH 2 



(21) 
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Notice that this time the dimension of the phase space 
has been increased from 3D to a 4-dimensional phase 
space spanned by the variables (x,y,z 7 v). However, as 
we will show quite soon, we can skip one of these variables 
and, consequently, the phase space corresponding to the 
present case can be reduced also to a 3-dimensional one. 

After the above choice of variables of the phase space, 
we can write the Friedmann constraint in terms of 
x,y,z,v: 



Sir. 



1 + x — z - 



vy-2 

v-1 



y/ v 2 y 2 — 2v(v — l)xy — 2vy ^ 
v — 1 

where now we have conveniently defined 



Slcurv : (^^) 



Sl r . 



6H 2 f R (l+4AH 2 



I J 6H 2 f R 



Notice that, contrary to the former case, here the effec- 
tive matter density parameter Sl m can be negative when- 
ever v < 1. In this case constraints on the domain of 
x,y,v arise from non- negativity of the expression under 
the square root: 



y> 



2(v-l)x + 2 



Other magnitudes of relevance for the analysis that 
can be straightforwardly put in terms of the above vari- 
ables are the effective equation of state (EOS) parameter 
uj e ff = — 1 — 2H/3H 2 , and the deceleration parameter 
q = -(l + H/H 2 ): 



1 - 2x 

ueff = —5—, q 



1 



respectively. 

Following the same procedure as before, after a bit of 
algebra, we can obtain the following autonomous system 
of ordinary differential equations out of the cosmological 
equations (fl"2" ]l , ([15 ]) : 



x' = -2x(x - 2) + 



2(x-2) 2 



v-1 

z[vy - 2(v - l)x - 2] 



y' = -2y(x-2) 



z' = z 2 — xz + 3z + 



v - 1 



2v{x - 2) 



v- 1 
2z(x-2) 
v — 1 



v' = -2v(x- 2) . 



(23) 



*s = {(%, V, z, v) : y > , v > , 

vy>2{v-l)x + 2}. (24) 

Worth noticing that in the variables x, y, z, v, the par- 
ticular case Q = is recovered from the present case in 
the formal limit v — > oo. 



2. Parametric ASODE 

It is noticeable the similarity between 2nd and 4th 
equations in (f2"3"|). However, by looking at the definitions 
of the variables y and v one sees that this is not casual. 
Actually, from (|2"Tj) one can see that 



y = f-iv , /x 



3a 2 ' 



(25) 



Hence, the ASODE (|23|) is not really a 4-dimensional 
system of equations but a 3-dimensional one. In fact 
it is a 3D, one-parameter ASODE. In the theory of the 
dynamical systems the parameter fi = 4/3 /3a 2 is called as 
"control parameter" . It usually controls the occurrence 
of bifurcations - a change of the value of the control 
parameter produces a change in the topology of the phase 
portrait. 5 That such bifurcations arise in the present 
study will be clear when we found that the asymptotic 
structure of the DBI-EH model (Q = => fi = 0) and 
that of the general, ghost-free case (fi ^ 0, fi > 0), are 
qualitatively different. 

In what follows we will skip the variable v, and will 
write the 4D autonomous system (|23l) as a 3D, one- 
parameter ASODE in the variables x, y, z: 



x = —2x(x — 2) 



2fx(x- 2) 2 

y - m 

z[y 2 - 2xy + 2fx(x - 1)] 



1/ — A* 



-2y(x-2) , 



z 1 = z 2 — xz + 3z 



2y(x - 2) 



y- n 

2nz(x-2) 
— — + 

y - m 



3Q>^(x,y,z), (26) 



where 



SV^(x,y,z) 



- X — Z h 

y\/y 2 - %xy + 2fi(x -~T) 



(27) 



The 4D phase space relevant for the cosmological 
model of interest in this case: 



5 Boundaries between regions of phase portrait with different 
topology are called as "bifurcation sets" . 



7 



TABLE III: Properties of the equilibrium points of the autonomous system (|26[) . We have defined the parameter y* — 
3/2 + fc/2 - (4/Lt- 9)/2fc, where fc = (2^ 2 - 18^ + 27 + 2^/^i 4 - 2/i 3 ) 1/3 . 
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-1 


-1 


ft 


2 


2/* 







fc > (ju > 2.25) 
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-1 


-1 



TABLE IV: Eigenvalues of the linearization matices corresponding to the critical points in table |TTTJ For the equilibrium point 
P-j the analytic expressions of the eigenvalues are extremely huge and complex, so that we have decided to show numeric results 
for several values of the control parameter fj, instead. 



Equilibrium Point fi Ai A2 A3 



ft 
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ft 
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-1 
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ft 
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-3/4 + iVfl/4 
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< fi < 2 


-2 


-4 


-4 


Pt 


< fi < 2 


4 


6 


-4 


ft 


5 


2.1 


-4.03 + 0.78i 


-4.03 - 0.78i 




10 


3.26 


-4.63 + O.Si 


-4.63 - 0.8i 




15 


4.19 


-5.1 + 0.55i 


-5.1 + 0.55i 



The one-parameter, 3D phase space spanned by the vari- 
ables x,y,z - relevant for the study of the asymptotic 
properties of the cosmological model (|12I13[) - can be 
finally defined in the following way: 



= {(x,y,z) :y>0, y 2 -2xy+2 A< (x-l) > 0} . (28) 
We want to emphasize that the bound 



B{x,y :f i) = y 2 -2xy + 2fi{x-l) >0, (29) 

has to be necessarily satisfied in order for the density 
parameter to be a real quantity. 



3. Equilibrium points 

There are found eight equilibrium points corresponding 
to the autonomous system of ODE (|26|) in the phase space 
Their most relevant properties are listed in the table 
IIIII while the eigenvalues of the linearization matrices 
corresponding to these points - invaluable as they are to 
judge about their stability -, are displayed in table H*Vl 

Notice that the point P3 in Tab. IIIII is not really a 
fixed point of (f2"FJl) , since the bound (|2"9"]) is not satisfied. 
Actually, for P 3 = (1/2,0,-3); 

= -fi < , 

since we are considering positive /i-s. 

Fixed points Pi, P2, and P4 are associated with in- 
finitely large values of the Hubble parameter H — > 00, 
while points P5-P7 correspond to de Sitter cosmological 
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FIG. 2: Drawing of the phase space to show the "multi- 
attractor" structure (ft = 1). This future asymptotic struc- 
ture is better seen in the r-fiow drawings in the lower panels. 



phases with (different) finite Hubble rates. 

The main asymptotic properties of the model with 
f(R) = 2(1 — y/1 — aR + j3Q) can be summarized as it 
follows: 

• The equilibrium point Pi is associated with 
curvature-dominated {Cl curv — 1, Q m = 0), non- 
inflationary dynamics as long as q = 0. It is the 
past attractor in Vf^. 

• Another non-inflationary phase can be associated 
with the fixed point P%, corresponding to matter 
domination (fi m = 1, Cl curv =0). It is a saddle 
equilibrium point in tyg. 

• The scaling solution fi m /fi curtI = —5/3 - point 
P3 - is associated with decelerated expansion (q = 
1/2). The effective "curvature" fluid mimics dust 
{w e ff = 0). However, as already noticed this is not 
a fixed point of ([26]) since it does not satisfy the 
bound P5]>. 

• Super-inflationary dynamics can be associated with 
the equilibrium point P4 in Tab IIIII (also a scaling 
solution since fl m /D, curv = —4/3). Actually, in this 
case w e ff = —5/3 < —1 (q = —2). Otherwise, the 
effective fluid mimics phantom behavior, which can 



be explained here as a curvature effect (no phan- 
tom fields at all). Point P4 is a future attractor 
in (this time not a global one, see below). It 
corresponds to a singular state (infinite curvature 
R — > 00) that could be identified with a catas- 
trophic fate inherent in a large class of phantom 
models. 

• Curvature-dominated (Cl curv — 1, fi m — 0) equi- 
librium points P^ 1 - corresponding to de Sitter 
inflationary phases of the cosmic expansion since 
q = — 1 - are also future attractors in These 
points exist whenever fi < 2. 

• Points P^ - associated also with de Sitter (infla- 
tionary) expansion - are saddle fixed points in the 
phase space (|2"8")l . These can be linked with scal- 
ing solutions since £l m /Q, curv = —4/5. The fact 
that the effective dimensionless matter energy den- 
sity parameter fi m is a negative parameter, is not 
a problem and is easily explained as due to our def- 
inition: 

fi = f — ^ Prn =( V \ Pm 

Actually, for y < /i (v < 1) fi m < 0. Points P^ 
exist only if fj, < 2. 

• It is found, additionally, another de Sitter sad- 
dle equilibrium configuration (point P7) in tyg. 
The fixed point P? can be associated also with a 
curvature-dominated solution (Cl curv = 1, fi m = 0) 
of the cosmological equations (fT2")l , ([TBI . This equi- 
librium point exists whenever the following inequal- 
ity holds: 

(2fi 2 - 18fx + 27 + 2vV - 2fi 3 ) 1/3 > 0. 

Numerical investigations show that this bound is 
equivalent to the following restriction on the values 
the control parameter /x can take: /1 > 2.25. 

We want to emphasize several interesting features of 
the asymptotic structure of the model (|12I13[) . Perhaps 
the most interesting fact is the "multi-attractor" struc- 
ture into the future of the r-development of ([2T)]) for val- 
ues of the control parameter /1 < 2. In this case there co- 
exist three future attractors (equilibrium points P4, P$~ , 
and -P 5 ~). Otherwise, for this interval of control parame- 
ter space there is no a unique global late-time attractor. 
Anyway, judging by the magnitudes of the eigenvalues, 
convergence of phase space trajectories towards the point 
P4 is stronger than towards points P 5 ± . 

Not of less importance is the existence of a saddle equi- 
librium point correlated with matter dominance - point 
P2 (^m = 1, Qcurv = 0). The existence of this point is 
crucial to explain the formation of structure. Existence 
of an inflationary de Sitter saddle equilibrium point P7 
could be important to explain early time inflation. 
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FIG. 3: Probe paths in phase space originated 
from different sets of initial conditions for f(R) = 
f 1 — \/T— aK 2 R + f3Q^\ . The control parameter 

H = 4/3 /3a 2 has been chosen to be fi — 5 so that points 
Pg and P^ do not exist. The flux in time is shown in the 
lower panels. The trajectories in converge towards the 
global ("super-inflationary") future attractor - equilibrium 
point P4 — (3, 0, 2) in Tab. Mil The past attractor is the 
curvature-dominated (non-inflationary) solution - point 
Pi = (1,0,0) in Tab. ITTT1 



Last but not least: the super-inflationary attractor 
point Pa (a global future attractor if fi > 2) represents 
another example of how to mimic catastrophic phantom 
behavior (w e ff < — 1) without any phantom fields. It is 
evident that in the model crossing of the phantom divide 
is possible (it is in fact unavoidable for fi > 2). 



4- Bifurcations 

To get additional information on the asymptotic dy- 
namical properties of the model under study one can 
rely on drawings of the phase space, where probe paths 
evolved from given (arbitrary) initial data probe the 
phase space and, eventually, uncover the nature of the 
relevant equilibrium points. Here, as in the former sub- 
section, we make a combined use of the standard tools 
of the (linear) dynamical systems analysis and of phase 
space drawings. This time the phase space drawings will 
serve as a graphic illustration of what it has been said 
in the former subsection (what is shown in tables Tab. 
IIIII and Tab. IIV[) . They clearly illustrate, in particular, 
bifurcations in respect to the control parameter /i. 

In the upper part of figure [5] probe trajectories origi- 
nated from different initial data, with control parameter 
fj, = 1, converge (increasing r) into three different at- 




FIG. 4: Plots of the phase space coordinate y vs fi for the 
equilibrium points P* and P^ (jj, < 2) - left panel. In the 
right panel it is shown a plot of the effective dimensionless 
matter density parameter Q. m vs /i for points P* (positive 
half), and P^ (negative half). At /i = 2 points Pg~ and P 5 ~ 
(also Pq and P 6 ~) coincide, so that there are 7 equilibrium 
points in tyg instead of 9. For fi > 2 these points do not exist 
at all (only 5 equilibrium points are found). This is a clear 
illustration of the bifurcation in the control parameter. 



tractor points: P 5 = (2,3.41, —1) (bunch of trajectories 
closer to the reader), P 5 ~ = (2,0.6,-1) (bunch to the 
center of the drawing), and P4 = (3,0,2) (set of trajec- 
tories at the end). This is more evident in the r-flow 
diagrams in the lower part of the figure. In the left panel 
the convergence towards three different points is appar- 
ent. In the right panel only two points seem to be points 
of convergence of the phase trajectories. Note, however, 
that in this diagram only the coordinate axes x and z 
are displayed. Hence, since the points P 5 + and P 5 ~ differ 
only in the coordinate y, in this drawing they appear as 
a single convergence point (the one closer to the reader) . 

In figure [3] a different value of the control parameter 
fi = 5 has been chosen. In this case, since /i > 2, equilib- 
rium points P 5 ± (also Pj^) do not occur. Consequently 
there is only a (global) future attractor: the point P4. As 
seen from Fig. [3] (upper part) probe paths in 'F^ converge 
into the past (decreasing r) towards the past attractor 
(equilibrium point Pi), while these converge into the fu- 
ture towards the attractor point P4. In the right panel 
a view from a different angle is shown. This asymptotic 
structure is confirmed by the r-flow diagrams in the lower 
part of the figure. 

The above figures (Fig. [2] and Fig. [3]) illustrate the 
bifurcation in respect to the control parameter /i. De- 
pending on the magnitude of the parameter fj, there are 
9 (/x < 2), 7 (;U = 2), or 5 (ji > 2) equilibrium points 
in the phase space 'Fg. Therefore, /x = 2 is the critical 
value or the bifurcation point. 

In the figure Fig. [4] we show bifurcation diagrams y 
vs /i - left panel, and O m vs ji — right panel, at the 
equilibrium points P^ 1 and P^. It is apparent that as 
fi — > 2 => y — > 2, points P 5 + and P 5 ~ tend to be a single 
equilibrium point (P5 = (2, 2, —1)). The same is true for 
the points P§ and P 6 ~ : at the critical value \x = 2 these 
are a single fixed point P^ = (2,2,3). In this case the 
four equilibrium points P^ 1 and P^ 1 reduce to just two 
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equilibrium points P5 and P§. For values fj, > 2 (not 
drawn in the figure) even P 5 and Pq are erased from the 
phase space. 

In the present model bifurcations in the control pa- 
rameter space demonstrate the strong dependence of the 
dynamic properties of the model on the values of the 
overall parameters a and /?. 



IV. DISCUSSION 

Knowledge of the equilibrium points in the phase space 
corresponding to a given cosmological model is a very 
important information since, independent on the initial 
conditions chosen, the orbits of the corresponding au- 
tonomous system of ODE will always evolve for some 
time in the neighborhood of these points. Besides, if the 
point were a stable attractor, independent on the initial 
conditions, the orbits will always be attracted towards 
it (either into the past or into the future). Going back 
to the original cosmological model, the existence of the 
equilibrium points can be correlated with generic cosmo- 
logical solutions that might really decide the fate and / or 
the origin of the cosmic evolution. In a sense the knowl- 
edge of the asymptotic properties of a given cosmological 
model is more relevant than the knowledge of a particular 
analytic solution of the corresponding cosmological equa- 
tions. While in the later case one might evolve the model 
from given initial data giving a concrete picture that can 
be tested against existing observational data, the knowl- 
edge of the asymptotic properties of the model gives the 
power to realize which will be the generic behavior of the 
model without solving the cosmological equations. In the 
dynamical systems language, for instance, a given partic- 
ular solution of the Einstein's equations is just a single 
point in the phase space. Hence, phase space orbits show 
the way the model drives the cosmological evolution from 
one particular solution into another one. Equilibrium 
points in the phase space will correspond to solutions of 
the cosmological (Einstein's) equations that, in a sense, 
are preferred by the model, i. e., are generic. The lack of 
equilibrium points that could be correlated with a given 
analytic solution of the model, amounts to say that this 
solution is not generic, otherwise it can be attained un- 
der a very carefully arrangement of the initial conditions 
only. 

The asymptotic structure of the model ([121131) (or 
(|8l9p ) is extremely rich and complex. Contrary to what 
is believed (see, for instance, Ref. [HI), there are found 
in this kind of f(R, 0)-theory, equilibrium points that 
can be associated with matter-dominated and matter- 
curvature scaling solutions. This is a very nice feature 
that makes the model attractive to do cosmology, since, 
by appropriately arranging the free parameters, there is 
room to accommodate the amount of structure obser- 
vations confirm it exists. The fact that, for /1 > 2.25 
(see existence conditions in Tab. IIII[) . two inflationary 
solutions with different Hubble rates, and with different 



stability properties, co-exist together: i) the inflationary 
de Sitter phase linked with the saddle fixed point P 7 ; 

jj2 _ 

3a(fc 2 + 3fc-4^ + 9) ' 

where k = (2^ 2 - 18/x + 27 + 2 vV - 2^ 3 ) 1 / 3 , and ii) the 
super-inflationary ghost-like phase (w e ff < — 1, q = —2) 
- point P4 - with infinite Hubble rate, is encouraging 
as well since, one can try to accommodate both early 
and late time inflationary stages into a united picture, 
where inflation is driven by curvature effects. In fact, the 
latter solution can be associated with late-time speedup 
of the cosmic expansion, while the former one might be 
linked with early inflation as long as, due to its stabil- 
ity properties, it is a transient phase. By appropriately 
choosing the value of the free parameter fi = 3a 2 /4:f3, 
one can manage to produce the necessary amount of in- 
flation at early times, besides, since P7 is a saddle critical 
point, then, exit from this inflationary stage is natural. 
Actually, fueled by the stability properties of the cor- 
responding cosmological phase, the dynamics forces the 
model to leave the saddle critical point Pj to, eventu- 
ally, approach to the global late-time super-inflationary 
attractor P4 (recall that we are considering fi > 2.25, 
so that the de Sitter attractors P^ do not exist). 6 For 
(i < 2, the possibility to explain early and late time in- 
flation in a united scheme is based on the co-existence of 
the super-inflationary attractor P4 and of the inflationary 
(de Sitter) saddle critical points P*- 7 

The existence of the non-inflationary, curvature dom- 
inated equilibrium point Pi = (1,0,0), being always the 
past attractor in the phase space, allows us to associate 
the starting point of every probe path in the phase space 
with a state characterized by infinitely large values of the 
Hubble parameter, a state usually linked with the big 
bang singularity. Therefore, contrary to what one should 
expect, the initial singularity is a generic feature of the 
model sourced by the action ([I]). 8 In case such a singu- 
larity were removed by the non-linear DBI dynamics, as 
expected, the starting point of the cosmic dynamics were 
to be linked with a regular solution resembling a finite 
big-bang event, which does not seem to be the case. It 
has to be pointed out, however, that, even if the generic 
behavior is associated with a big-bang singularity, there 



° Similar results about the possibility to have both inflationary 
quintessential and super-inflationary ghost-like phases in a model 
have been obtained for more elaborated models [2C|. 

7 Notice that, even if there co-exist de Sitter attractors - points 
Pg , and de Sitter solutions associated with saddle points P s , 
there is no room for a united description of early and late time 
inflationary stages of the cosmic expansion in this case, since the 
Hubble rate is the same for both asymptotic stages. 

8 This kind of statements has to be taken with caution since the 
dynamical systems analysis uncovers the asymptotic properties 
of the dynamics. Hence, for instance, the asymptotic state where 
H 2 blows up might be reach in an infinite time, in which case one 
would not want to say that the space-time is properly singular. 
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might exist particular solutions that are regular at the 
beginning of the cosmic expansion. To show that the lat- 
ter statement is indeed true, let us write the combination 
a(R + AQ) in the action (|2J[3|), in terms of the variables 
of the phase space: 

a(R + AG) — = ■ 

There is a class of trajectories in the x — y phase plane: 

m . y 2 f(y) 

x = 1 , 

where f(y) is a generic function such that, in the limit 
y — > 0, it is finite (note from the above expression that, 
at y = 0, x = 1). For these particular trajectories as one 
approaches the point Pi — (1,0,0), the above curvature 
invariant 

a(R + AG) = -f(y), 

so that, at x = 1, y = it is finite, otherwise, the action 
is finite and there is no singularity associated with the 
past attractor P\. 

V. CONCLUSION 

In the present paper we have performed a through- 
out study of the asymptotic properties of a cosmological 
model that is based on a non-linear (DBI-type) modi- 
fication of Einstein-Hilbert general relativity previously 
proposed in |9|. The modification includes, in addition 
to the standard Einstein-Hilbert term within the square 
root in the action, also the Gauss-Bonnet invariant. The- 
ories of this kind are oriented to remove the pathological 
spacetime singularities from the original EH theory. In 
former works in the bibliography the present model has 
been demonstrated to be ghost-free and, for a given in- 
terval of the space of parameters, it is also free of the 
Coulomb-like Schwarszchild sing ularity @. 

The study has been performed in two steps: i) first we 
explored the asymptotic properties of a simpler model 
where the Gauss-Bonnet invariant is absent, and then ii) 
we studied properly the model of interest. The idea was 
to solve first a related case easier to handle so that we 



could get some insight into the more complicated situa- 
tion. Even if the DBI model without the GB invariant 
has a ghost (spin-zero) degree of freedom |l8|, so that it 
is not field-theoretically consistent (negative-norm states 
arise when invoking quantum computations), the study 
showed, nonetheless, that even this simpler - theoreti- 
cally inconsistent - model has very nice asymptotic fea- 
tures that make it attractive to do cosmology. 

In general the asymptotic structure of the Friedmann- 
Robertson- Walker cosmological model originated from 
the action (TTJ) is extremely complex and rich. There 
can be found equilibrium points associated with matter- 
dominated, matter-curvature scaling, de Sitter and even 
phantom-like solutions. The multi-attractor structure in 
the future r-development of the autonomous system of 
ODE (|26p is a very interesting finding. Even more in- 
teresting are the bifurcations in the control parameter p,- 
space, showing the strong dependence of the model upon 
the values of the overall parameters a, and /3 (this is very 
convenient as long as makes the model very flexible to ac- 
commodate different kinds of cosmological dynamics). 

Our results concerning the occurrence of a future sin- 
gularity as well as of a past big-bang singularity, confirm 
the results of previous studies on this subject within a 
large class of modified theories of gravity [Hj . Otherwise 
stated, contrary to what one should expect from a non- 
linear modification of general relativity such as the one 
based on the action (fTJ) , removal of the big bang singular- 
ity is not a generic feature of the corresponding cosmo- 
logical model. Instead, the starting point of the cosmic 
expansion is a state with infinitely large values of the 
Hubble parameter squared H 2 — > oo, which is associated 
with a past attractor in the phase space. Additionally, 
there is a future super-inflationary attractor associated 
with a finite-time singularity of the kind inherent in a 
large class of phantom models of dark energy. 
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